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Abstract. X-ray observations of middle-aged pulsars allow one to study nonther- 
mal radiation from pulsar magnetospheres and thermal radiation from neutron 
star (NS) surfaces. In particular, from the analysis of thermal radiation one can 
infer the surface temperatures and radii of NSs, which is important for investi- 
gating evolution of these objects and constraining the equation of state of the 
superdense matter in the NS interiors. Here we present results of XMM obser- 
vations of three middle-aged pulsars, J0538+2817, B0656+14 and J0633+1746 
(Geminga), and briefly discuss mechanisms of their X-ray emission. 
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1. Introduction 

A large number of isolated neutron stars 
(NSs) of various types, including active 
radio pulsars, have been observed with 
the currently operating large X-ray ob- 
servatories, Chandra and XMM. Recent 
reviews on Chandra observations of NSs 
can be found in Pavlov et al. (2002; 
[PZS02]) and Weisskopf (2002), while first 
XMM results are described by Becker & 
Aschenbach (2002). Here we discuss XMM 
observations of three pulsars, J0538+2817, 
B0656+14 (J0538 and B0656, hereafter) 
and J0633+1746 (Geminga), with ages of 
30-300 kyr and spin-down energy loss rates 
of (3-5) x 10 34 erg s^ 1 (see Table 1). 

X-ray emission from an active pulsar 
generally consists of two components, ther- 
mal and nonthermal. In very young, bright 
pulsars (r < 10 kyr) the thermal compo- 



nent is buried under the powerful nonther- 
mal emission, while X-ray emission from 
old pulsars (r > 10 3 kyr) is usually very 
faint. The intermediate group of middle- 
aged pulsars is particularly interesting for 
X-ray studies because they are expected 
to be bright and exhibit both thermal and 
nonthermal components. Below we present 
results of spectral and timing analysis of 
three middle-aged pulsars observed with 
the EPIC-pn instrument on board XMM. 



2. PSR J0538+2817 

This radio pulsar was observed with XMM 
in March 2002 (11.3 ks effective exposure 
for EPIC-pn, in small- window mode). The 
detected spectrum of J0538 can be fitted 
with a single blackbody model of T££ ~ 
2.12 MK, without invoking a nonthermal 
(power-law) component (McGowen et al. 
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ms 


kyr 


10 34 ergs -1 


10 12 G 


kpc 


10~ 12 cgs 


10 32 ergs^ 1 


10 32 crgs^ 1 


J0538 


143 


30 
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1.2 


0.6 


6.5 


11.2 


B0656 


385 


110 


4 


5 
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7.1 


1.3 


1.9 


Geminga 


237 


340 


3 


2 


0.2 


1.8 


0.2 


0.5 



Table 1. Pulsar parameters: period P, age r (P/2P for B0656 and Geminga, and true 
age for J0538 — Kramer et al. 2003), rotational energy loss E, surface magnetic field B, 
and distance d. Also, the energy flux Fx in 0.2-10 keV (as detected with EPIC-pn), the 
corresponding luminosity Lx (corrected for interstellar absorption) and the bolometric 
thermal luminosity Lboi are given. 



2003). The blackbody model gives a small 
X-ray emitting area, with a radius ~ 
1.7 km 1 , much smaller than the expected 
NS radius ~10 km. Such an area might be 
interpreted as a hot spot on the pulsar's 
surface. An alternative interpretation in- 
volves NS atmosphere models (see Zavlin & 
Pavlov 2002 [ZP02] for a review). A hydro- 
gen atmosphere model with B = 1 x 10 12 
G fits the data even better (Fig. 1) and 
yields a lower (unredshifted) surface tem- 
perature T c ff ~ 1.12 MK and a radius 
R ~ 10.5 km, consistent with the expected 
NS radius. An upper limit on the non- 
thermal X-ray luminosity in 0.1-10 keV is 

jnonth < 1 x 1Q 31 erg s -l = 2 x IQ-^E. 

Similar results were obtained by Romani & 
Ng (2003) from a Chandra observation, al- 
beit with somewhat different model param- 
eters (the discrepancy is probably due to a 
pile-up effect in the Chandra data). They 
also found an indication of a pulsar-wind 
nebula around J0538, with a size of ~ 18", 
that cannot be resolved in the XMM data. 
However, the estimated luminosity of the 
nebula, ~ 6 x 10 _5 .E, is too low to seriously 
affect the analysis of the EPIC-pn data. 

The EPIC-pn data also show pulsations 
of the X-ray flux, with one broad, asym- 
metric pulse per period (pulsed fraction 
f p ~ 25%; McGowen et al. 2003). The 
phases of pulse maxima at energies below 
and above 0.8 keV differ by - 75° (Fig. 
2). Such pulsations show that the thermal 



1 All radii and luminosities are given for the 
distances listed in Table 1. 
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Fig. 1. EPIC-pn count rate spectrum of 
J0538 fitted with a magnetized hydrogen 
atmosphere model. 

emission is intrinsically anisotropic and in- 
dicate strong nonuniformity of the surface 
temperature and magnetic field. 

3. PSR B0656+14 

B0656 was observed with XMM in October 
2001, with 6.0 ks and 25.0 ks EPIC-pn ex- 
posures in small- window and timing modes, 
respectively. The results of the spectral 
analysis of the EPIC-pn data are quite con- 
sistent with those derived from Chandra 
observations (PZS02). The X-ray spectrum 
shows at least two distinct components, 
thermal and nonthermal. The best fit is 
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2. X-ray light curves of J0538. 
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Fig. 3. EPIC-pn count rate spectrum of 
B0656 taken in two observational modes 
and fitted with a three-component model. 



provided by a model consisting of soft (TS) 
and hard (TH) blackbody components, and 
a power law (PL; Fig. 3). The model pa- 
rameters are T b °g s ~ 0.82 MK, R™ „ ^ 7.3 



km (for TS), TJ 



bb.h 



L bb,s 

1.72 MK, 



0.5 km (for TH), and the photon index 
7 ~ 1.5 (for PL). This model fits well 
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Fig. 4. Multiwavelength spectrum of 
B0656, with extrapolations of the total 
X-ray spectrum (long dashes) and its three 
components into the optical domain. 



the multiwavelength spectrum of B0656, 
from IR through X-ray energies (Fig. 4; see 
also PZS02 for details). However, even the 
larger radius, R = [1 + z]~ 1 R^ ) S w 5-6 km 
(z is the gravitational redshift), is too small 
to be regarded as a NS radius. Using NS 
hydrogen atmosphere models for the ther- 
mal component yields distances d ~ 0.1 kpc 
(at R = 10 km), a factor of three smaller 
than actually measured. The nonthermal 
component contributes only 1% of the to- 
tal luminosity in 0.1-10 keV, L^ onth ~2x 
10 30 crgs -i = 5 x iq-^e. 

The timing analysis of the EPIC-pn 
data confirms that both the pulse shape 
and pulsed fraction strongly depend on 
photon energy (Fig. 5). The X-ray pulse 
profile of B0656 apparently consists of two 
components (centered at cf) ~ 0.6 and 1.0), 
whose relative contributions depend on en- 
ergy. The nonthermal emission at E > 2 
keV shows one narrow peak per period, 
with f p > 30%. Unfortunately, strong back- 
ground contamination (especially in timing 
mode) prevents a more accurate analysis of 
pulsations at these energies. 
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Fig. 5. X-ray light curves of B0656. 



4. Geminga 

A long XMM observation of this famous 
7-ray pulsar in April 2002 revealed a bow- 
shock nebula caused probably by the su- 
personic motion of the NS (Caraveo et 
al. 2003). The EPIC-pn effective exposure 
(small-window mode) was 71.4 ks. A large 
number of photons collected allows one to 
accurately measure the spectral parameters 
of the pulsar. The phase-integrated spec- 
trum nicely fits a blackbody model, with 
T b °g ~ 0.50 MK and R^ h ~ 10.6 km, plus a 
power-law component with 7 ~ 1.85 (Fig. 
6). Similar to B0656, the NS atmosphere 
models yield too large radius-to-distance 
ratios. The nonthermal component pro- 
vides 7% of the total luminosity in the 0.1- 
10 keV range 2 , i^ onth ~ 3 x 10 30 < 



ergs 



2 The nebula might contribute a small frac- 
tion in the estimated nonthermal flux of the 
pulsar; a quantitative estimate is expected 
from a Chandra observation. 
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Fig. 6. EPIC-pn count rate spectrum of 
Geminga fitted with a two-component, 
blackbody (BB) plus power-law (PL), 
model. 
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Fig. 7. Multiwavclcngth spectrum of 
Geminga, with extrapolations of the to- 
tal X-ray spectrum (long dashes) and 
its thermal component into the optical 
domain. 



1 x 10~ 4 i?. The composite X-ray spectrum 
extrapolated into the optical domain (Fig. 
7) exceeds the observed spectral fluxes by 
a factor of 50-100, whereas the extrapo- 
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Fig. 8. X-ray light curves of Gcminga. 

lated thermal component is close to the 
Rayleigh- Jeans spectrum observed in UV 
(Kargaltsev et al. 2004). 

Similar to the case of B0656, the X-ray 
pulse profile shows a complex energy de- 
pendence (Fig. 8), with f p growing up to 
almost 100% at E > 2 keV, where the emis- 
sion is purely nonthermal. The shape of the 
light curve extracted in the 2-4 keV range 
is similar to that observed in gamma-rays, 
at E > 100 MeV (Mayer-Hasselwander et 
al. 1994). The two narrow peaks of the 
nonthermal emission are superimposed on 
the broad peak at lower energies. Phase- 
resolved spectroscopy shows that the pho- 
ton index varies with phase from 7 = 1.7 
at <f> ~ 0.4 (near the minimum of the com- 
posite light curve) to 7 = 2.0 at the phase 
<p ~ 0.7 of the broader peak in the 2-4 keV 
range. 

5. Discussion 

The main contribution to the X-ray fluxes 
of the three pulsars is provided by ther- 
mal emission from the their surfaces. The 
younger age of J0538 naturally explains 
its higher surface temperature, that is 
consistent with the standard NS cooling 



scenario. As discussed in ZP02, at such 
high temperatures ( > 1 MK) and relatively 
low magnetic fields, the NS surface is ex- 
pected to be covered with a gaseous at- 
mosphere, strongly ionized if comprised of 
hydrogen. Contrary to J0538, the older 
B0656 and Gcminga are significantly colder 
and have stronger magnetic fileds. At such 
conditions, hydrogen atmospheres become 
denser and less ionized, and the currently 
available atmosphere models become less 
reliable. On the other hand, diffusive nu- 
clear burning can reduce the hydrogen 
abundance in envelopes of cooling NSs 
(Chang & Bildsten 2003). In this case, the 
cooling hydrogen-depleted envelopes may 
undergo a phase transition, condensing into 
a liquid or solid surface. This may explain 
why the fits with hydrogen atmosphere 
models do not provide reasonable param- 
eters for B0656 and Geminga, contrary to 
the younger and hotter J0538. 

The observed temporal behavior of X- 
ray fluxes from these three pulsars indi- 
cates that their thermal radiation is locally 
anisotropic, in obvious contradiction with 
the simplistic blackbody interpretation of 
the phase-integrated spectra of B0656 and 
Geminga. Moreover, the asymmetry of the 
soft X-ray pulses (particularly for J0538 
and B0656) hints that the surface distri- 
butions of temperature and magnetic field 
are not azimuthally symmetric, suggesting 
a strong multipolar component of the mag- 
netic field or a decentered magnetic dipolc. 
Whatever is the true nature of the com- 
plicated light curves, the NS parameters 
inferred from the blackbody spectral fits 
should be taken with caution. In particular, 
the comparison of the blackbody temper- 
atures with the predictions of NS cooling 
theories may not be reliable, and the black- 
body radii may be quite different from the 
actual NS radii. 

The X-ray spectra of nonthermal emis- 
sion from B0656 and Geminga are close to 
power laws, with a phase-dependent photon 
index in the latter case. While the power- 
law slopes and intensities (but not the light 
curves — Kern et al. 2003) in the opti- 
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cal and X-ray bands are in agreement with 
each other for B0656, they strongly dis- 
agree for Geminga. It might be explained 
by the presence of two different populations 
of relativistic electrons in the pulsar's mag- 
netospherc. The younger and more ener- 
getic J0538 shows no nonthermal X-ray ra- 
diation, in contradiction with a general ex- 
pectation that younger pulsars with larger 
E are stronger nonthermal emitters. A sim- 
plest explanation of this contradiction is 
that the X-ray pulsar beam cannot be seen 
from the Earth. Alternatively J0538 could 
be one of several pulsars with undcrlumi- 
nous magnetospheric X-ray emission, — for 
instance, the famous Vela pulsar, whose age 
and thermal emission are similar to those of 
J0538, shows L n ° nth /E - 3 x 1(T 6 (Pavlov 
et al. 2001), well below the upper limit of 
2 x 10~ 4 inferred for J0538 from the XMM 
data. 

Generally, the new observations confirm 
that middle-aged pulsars share a number 
of common features (in particular, pulsed 
thermal radiation), but their X-ray and, es- 
pecially, multiwavelength emission proper- 
ties are rather diverse. Thanks to the much 
improved quality of X-ray data, we are now 
quite certain that we do see the NS sur- 
face layers in soft X-rays, with tempera- 
tures < 1 MK, generally decreasing with 
NS age, and we do see strongly pulsed 
magnetospheric radiation at higher X-ray 
energies. On the other hand, we still do 
not completely understand the mechanisms 
of thermal emission (hence, cannot accu- 
rately measure the NS temperatures and 
radii, despite very small statistical errors), 
and we lack quantitative models for mag- 
netospheric emission (hence, cannot infer 
the complex physics of NS magnetospheres 
from the measured phenomenological pa- 
rameters). However, it is not surprising 
that we have not yet been able to solve 
all the puzzles from observations of a tiny 
fraction of these exotic (and diverse!) ob- 
jects. The outstanding capabilities of the 
currently operating observatories, in both 
X-rays and optical, allow observations of a 
much larger sample of NSs, including ther- 



mally emitting midllc-aged pulsars, and we 
expect that our understanding of these elu- 
sive objects will be much better in a few 
years, provided enough observational time 
is allocated to NS studies. 
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